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Abstract 

In this paper (the first of a series) we describe the construction of fixed 
point actions for lattice SU(3) pure gauge theory. Fixed point actions have 
scale invariant instanton solutions and the spectrum of their quadratic 
part is exact (they are classical perfect actions). We argue that the fixed 
point action is even 1-loop quantum perfect, i.e. in its physical predictions 
there are no g 2 a n cut-off effects for any n. We discuss the construction of 
fixed point operators and present examples. The lowest order qq potential 
V(r) obtained from the fixed point Polyakov loop correlator is free of any 
cut-off effects which go to zero as an inverse power of the distance r. 
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1 Introduction and summary 



Using a discrete space-time lattice as an ultraviolet regulator for a quantum 
field theory also opens the way to many non-perturbative methods. At the same 
time however the lattice introduces various artifacts like rotational symmetry 
breaking. Cut-off independent continuum quantities are obtained in the limit 
of zero lattice spacing. Depending on the relevance of the lattice artifacts this 
limit can be reached on relatively coarse lattices or one might need lattices where 
the lattice spacing in physical units is extremely small. The main problem for 
numerical lattice calculations is the control of the finite lattice spacing effects. 

The relevance of the lattice artifacts on the numerical results can strongly 
depend on the specific regularization or lattice action. Most numerical calcu- 
lations to date of QCD use the Wilson plaquette action. This action is the 
simplest one which can be used but it has some undesirable properties. For 
example, asymptotic scaling sets in slowly and non-smoothly. Even more im- 
portant is that the continuum limit is reached at lattice spacing less then 0.1 
fm and, consequently, nowadays typical pure gauge or quenched calculations are 
done on spatial volumes 32 3 or larger . 

Two of us Q suggested a radical solution, to use a perfect lattice action 
which is completely free of lattice artifacts. That such a perfect action exists 
follows from Wilson's renormalization group theory The fixed point (FP) of 
a renormalization group (RG) transformation and the renormalized trajectory 
(RT) along which the action runs under repeated RG transformations form 
a perfect quantum action. The FP action itself reproduces all the important 
properties of the continuum classical action (it is a classical perfect action). 
Although at finite coupling values the FP action is not "quantum perfect" ^ it 
is expected to be a very good first approximation. This program was successfully 
carried out for the 2-dimensional non-linear a model. An action parametrized 
by about 20 parameters was found that showed no finite lattice spacing effects 
even at correlation length £ ~ 3 lattice spacings. The computational overhead 
to simulate this action (about a factor of three relative to the nearest neighbor 
action) is negligible compared to the gain obtained by the great reduction of 
lattice artifacts. 

In this and a follow up paper g we report on a similar project for the 
4-dimensional pure gauge SU(3) theory. We construct FP actions and FP oper- 
ators and investigate their properties. In this first paper the basic equations are 
presented and the properties of the FP action and operators are investigated 
on fields with small fluctuations. The methods used are mostly analytic. In 

3 In some limiting cases the FP action might be trivially related to the renormalized tra- 
jectory |y. 
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Figure 1: The quantity rV(r) constructed from the correlators of the FP 
Polyakov loops and using the FP action. Note that I/At: = .0796. 

the second paper we turn to the problem of using the FP action in numerical 
simulations on coarse lattices. A simple parametrization is derived and the scal- 
ing properties of the action are studied by measuring torelon masses and the 
potential. 

In the present paper we consider two different type of RG transformations. 
Both transformations contain free parameters which are used to make the FP 
action as short ranged as possible. For any practical application this is a basic 
requirement. The FP action is defined by a classical saddle point equation. It 
follows from this equation that the FP action has exact stable instantons. The 
value of the action for the instantons is 87r 2 independently of the size of the 
instanton (scale invariance). The FP equation can be solved analytically on 
fields with small fluctuations by expanding in powers of the vector potentials. 
The quadratic part of the FP action describes massless gluons with an exact 
relativistic spectrum: E = \p\, pi G (—00,00). This is to be contrasted with 
that of the Wilson action where this dispersion relation is valid only for |g| Cl 
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Figure 2: The quantity rV{r) constructed for the Symanzik action. 



and, in addition, the spectrum is restricted to the Brillouin zone \qi\ < 7T. For 
small \q\ the Symanzik tree level improved action performs better, but 
then, in the middle of the Brillouin zone it breaks down completely producing 
complex energies. 

An exciting question is how the FP action performs in 1-loop perturbation 
theory. We present a formal RG argument that the physical predictions are free 
of cut-off effects 0(a n ) and 0(g 2 a n ) for any n provided the size of the system 
is larger then the range of the action. Is therefore the classical perfect action 
automatically 1-loop quantum perfect? Our RG arguments are formal, not 
rigorous. The statement is supported further by a detailed 1-loop calculation 
in the non-linear a model || where the cut-off effects of the mass gap are 
studied in a finite box of size L. The power like decaying cut-off effects (a/L) n 
were found several orders of magnitude smaller (and consistent with zero within 
the errors of the calculation) than those in the standard or in the tree level 
improved Symanzik model. 
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Figure 3: The quantity rV(r) constructed for the Wilson action. 



We need not only perfect actions but perfect Polyakov, or Wilson loops, 
currents, etc. as well. As a first step we discuss here how to construct FP 
operators. We derive the explicit form of the FP Polyakov loop on fields with 
small fluctuations and obtain the corresponding tree level qq potential V(r). We 
show that V(r) is free of any cut-off effects which decay as a power r~ k for large 
r. This property is to be contrasted with that of the standard and Symanzik 
action. In the latter case only the leading 0(et 2 ) distortion is canceled. The 
situation is illustrated in Figs. |l|, ^ and [| 

Kerres, Mack, and Palma|| have recently described the construction of a 
perfect action for scalar electrodynamics in four dimensions and at finite tem- 
perature. Some of the basic ideas of || overlap with those introduced in Q and 
in the present work. 

The outline of the paper is as follows. In Section 2 the two different RG 
transformations we considered are defined. In Section 3 the FP equation is 
derived and the statement about classical solutions of the FP action is discussed 
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briefly. In Sections 4 and 5 the FP action is constructed at the quadratic level 
and the spectrum is studied. In Section 6 the main steps leading to the cubic 
terms of the action are presented. Section 7 deals with the question of whether 
the FP action is 1-loop quantum perfect. The problem of FP operators is 
discussed in Section 8. 



2 The RG transformation 



Consider an SU(N) gauge theory defined on a hypercubic lattice. The partition 
function is 

Z = J ' DUe~^ u \ (1) 

where (3S(U) is some lattice representation of the continuum action. It is a func- 
tion of the products of link variables U^ L {n) — e %A ^ n ^ gSU(N) along arbitrary 
closed loops. The normalization is fixed such that on smooth configurations the 
action takes the standard continuum form 

0S(U) -» ^ / d4xTr ( F ^ F ^) > ( 2 ) 

where 

2N 

IF' (3) 



Let us denote the couplings of /3S(U) by /?, 02,03, . . .. The action can be 
represented by a point in this infinite dimensional coupling constant space as 
shown in fig. |J. 

Under repeated real space renormalization group transformations the action 
moves in this multidimensional parameter space. We shall consider RG trans- 
formations with a scale factor of 2. The blocked link variable V^(ns) SSU(N), 
which lives on the coarse lattice with lattice unit a' = 2a, is coupled to a local 
average of the original link variables. The new action is defined as 

e -fi's'{V) = J DU exp {-(3 (S(U) + T(U, V))} , (4) 
where the blocking kernel T(U, V) is taken in the form 

T(U,V) = J2 {-^RcTr(V/> s )Qt> B )) +M(Q^(n B ))} . (5) 



G 




g=l/p 



Figure 4: Schematic flow diagram for asymptotically free gauge theories. 

The N x N complex matrix (tib ) is an average of the fine link variables 
in the neighbourhood of the coarse link Is = (% , \x) . Its explicit form will 
be specified later. The last term in Eq. (Q) assures that the partition function 
remains invariant under the RG transformation, 

e 0M(Q) = J dWexp | J g^.Be'&(WQ t )} , W £ SU(N) . (6) 

The parameter k is free and will be used for optimization as in We require 
the averaging procedure to be gauge invariant 

T(U 9 ,V 9 ) =T(U,V), (7) 

where U 9 and V 9 are gauge transformed configurations on the fine and coarse 
lattice, respectively. This property, together with the gauge symmetry of the 
action S(U), assures that S'(V) is also gauge invariant. 

In this paper we consider two different block transformations which will be 
referred to as type I and type II, respectively. The first transformation we used 
is a modified version of the Swendsen transformation JTo| | shown in fig. ^| where 
c, the relative weight of the staples vs. the central link, is a tunable parameter: 
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' i ~) i . = (l"6c) o > o > o + C i o + ... 
Figure 5: The renormalization group transformation of type I. 

Q„{n B ) = (1 - 60)^(71)^(71 + im) + 
cJ2 [U v (n)Up(n + v)U ll (n + iM + v)Ut(n + 2iM)+ (8) 

Ul{n - v)U^{n - v)U^{n + fi - ^^(n + 2^ - 1/)] . 

The second transformation uses smearing to create fuzzy links. The ma- 
trix Qu(jib) is the product of the two fuzzy link variables obtained after two 
smearing steps. The construction is shown in fig. 0. 



(l-6d) 



(l-6d) 



Figure 6: The renormalization group transformation of type II. 

Here c and d are tunable parameters. With c = 1/7 the transformation 
Eq. (|J) is identical to the one used in the /3 function calculations earlier |l(]] . 



It is expected that these RG transformations have a fixed point (FP) on the 
(3 = 00 critical surface. The FP has one marginally relevant and infinitely many 
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irrelevant directions. The trajectory which leaves the FP along the marginal 
direction is the renormalized trajectory (RT). 

As it was argued in 0] the FP and the RT describe perfect actions. There are 
no lattice artifacts along the RT. For reasons discussed in the Introduction and 
later, we shall refer to the FP action (3S FP (U) as the classical perfect action. 
This is the action we shall consider in this work. Although this action does not 
lie on the RT at finite (3 values, it has amazingly nice properties at large (3 and 
shows no, or tiny, cut-off effects even on very coarse configurations. 



3 The fixed point equation 



On the critical surface (i — > oo Eq. (Q) reduces to a saddle point problem giving 
S FP (V)=mm(S FP (U)+T(U 1 V)), (9) 

where S FP is the FP action. For a given coarse configuration {V} there are 
infinitely many minimizing {U} configurations which differ by gauge transfor- 
mations leaving the averages Q^ub) invariant. The value of the FP action 
S FP (V) is not influenced by this freedom, however. 

Using Eqs. the saddle point equation can be written as 

S FP (V)=miJs FP (U)-^Y, [ReTr(V> fl )Q> B ))-/(Q> s ))]l, 

(10) 

where 

f(Q) = max{ReTr(VFQ t )} , W € SU(iV). (II) 



This equation has a very similar structure to the one in the d = 2 non 
linear cr-model and its implications for the classical solutions of S FP are 
the same: if the configuration {V} satisfies the FP classical equations and is a 
local minimum of S FP (V) then the configuration {U(V)} on the finer lattice 
minimizing the r.h.s.. of Eq. ( |Io| ) satisfies the FP equations as well. In addition, 
the value of the action remains unchanged, S FP (V) = S FP (U(V)). 

Since the proof goes along the same lines as in the d = 2 cr-model , we 
can be brief here. If {V} is a solution of the classical equations of motion 
5S FP /5V = 0, then Eq. ( |To| ) implies that on the minimizing configuration 
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{U(V)} the expression in square brackets in Eq. ( |10| ) takes its maximum value, 
i.e. zero. Then it follows that 

SS FP (U) 
SU 

and S FP (V) = S FP (U(V)). 

According to this result the FP action has exact scale invariant instanton 
solutions with action = 8tt 2 as in the continuum theory. This property allows 
one to define a correct topological charge on the lattice Jll]]. In this paper, 
however, we will not follow this problem further. 



= 0, 



(12) 



4 The FP action at the quadratic level 

The FP equation Eq. ([!(]) is valid for any configuration {V}, be it smooth or 
rough. It is a highly non-trivial equation. On {V} configurations with small 
fluctuations, however Eq. ( |To| ) can be expanded in powers of the vector potentials 
and studied analytically. 

On the quadratic level we can write^. 

2NS(U) = P^(r)1r (A„(n + r)A v (n)) + O (A 3 ) (13) 

n,r 

= ^ E M*)^ (M-k)A„{kj) + o (i 3 ) , 

where p^ v {r) (in momentum space p M „(fc)) are the quadratic couplings to be 
determinecj^J. Gauge symmetry and the normalization condition Eq. (Q) give 
the constraints 

fc>„(fc) = 0, (14) 
Pnu{k) — > k 2 5 MV - k^k v for k — > 0, 

where we introduced the notation = e lkfi — 1. For both RG transformations 
(fig. H and fig. ||) Q^{ub) is a sum over products of U matrices along paths 
connecting the points n and n + 2/z. It is easy to show that 

ReTr (^(n B )Qj> fl )) - /(Q„(n fl )) = -^Tr (r„(n fl ) - B„(n B )) 2 + O (cubic) , 
" (15) 

4 We suppress the index "FP". 

5 In the following we suppress the tilde for Fourier transformed quantities. 
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where 



r^(ns) = uj^(2n B - r)A v (r) . 



(16) 



In Eq. (plf) the tensor uj^ is fixed by the form of the RG transformation kernel 
in fig. gpand fig. ^: 



v (k) = (l + e lfc ") [S^ (l-c(kk*)) + 



type I, 



(17) 



w„„(fc) = (1 + e 1 ^) 



(l - d(kk*j) + dfo- d(kk*j) %k* v 



At the quadratic level the recursion relation has the form 
^Y,PU k ^(M-kB)BAkB)) = 



V B 



type II. 



(18) 



^\^Y,P^ k ^(M-k)Au(k)) + (19) 

I k 

E Tr K^C"**) - M-ks)) (T^{k B ) - B^(k B ))} 1 , 

where V B = j^V is the volume of the blocked lattice. To proceed further we 
have to introduce a temporary gauge fixing to get rid of the remaining gauge 
freedom at the sites of the fine lattice not representing a site on the coarse 
lattice. One obtains then from Eq. (JL9|) a recursion relation for the inverse of p: 



KM = ^ £ [ W (f + nl)D(f + ^ f (f + ttZ)] ^ + ~J,u , (20) 



1 = 



where / = (Iq, l±, li, I3) is an integer vector, the summation goes over ^ = 0,1 
and 

D^(q B ) = (p-\q B ))^ ■ (21) 



A way to find the FP solution of Eq. ( p0|) is to start from a D$ (q) with 
some gauge fixing parameter a 



(qq*) 



(qq* 



(22) 



and iterate Eq. (|2C|) to the FP solution . The inverse of has a smooth 
a — > 00 limit which corresponds to removing the gauge fixing. The resulting 
Puu(k) i s transverse and satisfies Eq. (fl4|). 
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r 
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1 


-0.6718 





-0.8007 


2 
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1 


-0.03100 
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-0.00422 


2 


-0.00170 


10 


-0.05951 


10 


0.01129 


11 


-0.04281 


10 1 


-0.00208 


10 1 


0.01614 


10 11 


0.00880 


111 


-0.02724 


2 10 


-0.00387 



Table 1: Some of the elements of poo(^) and pio(r) for the RG transformation 
of type I. 

The program sketched above is not completely trivial. The main steps and 
some intermediate results are collected in Appendix A. 

The fixed point p^f depends on the parameters of the RG transformation 
(k, c), or (k, d) for the kernels of type I and type II, respectively. Although 
none of the theoretical properties of the FP action depend on these parameters, 
the range of p^f (r) (the range of the interaction) changes significantly as these 
parameters are varied. It is vital in applications that the interaction is short 
ranged. Using this requirement we found that the following parameter values 
are close to optimal: 

K = 12.Q, c = 0.12, type I, (23) 
k=10.5, d= 0.077, type II. 

Using these values the significantly non-zero couplings lie in the unit hypercube 
and puv( r ) S oes rapidly to zero with the distance. The decrease seems to be as 
fast as ~ exp(— 3r). In Table [j] and || we collected some of the elements of p^„(r) 
for the RG transformation of type I and II, respectively. Using these tables, plus 
the cyclic permutation of the coordinates and the reflection properties of p My (r) 
(Appendix B), further elements can be obtained. 



5 The spectrum of the quadratic FP action 

Although the quadratic FP action Eq. ( ]l3] ) lives on the lattice, it describes 
transverse massless gluons with an exact relativistic spectrum E{k) = \k\, 
hi 6 (—00,00). This follows from the fact that the FP action represents in 
fact an infinitely fine lattice, i.e. the continuum. A similar property has been 
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r 


Poo(r) 


r 


Pio{r) 


1 


-0.5872 





-0.7777 


2 


0.00769 


1 


-0.04813 


3 


-0.00035 


2 


-0.00081 


10 


0.00467 


10 


-0.00102 


11 


-0.07964 


10 1 


-0.00063 


10 1 


-0.00298 


10 11 


0.00804 


111 


-0.02490 


2 10 


-0.00571 



Table 2: Some of the elements of poo(^) and pio(r) for the RG transformation 
of type II. 

demonstrated earlier for the d = 2 non-linear a-model where this state- 
ment follows trivially from an explicit closed representation of the FP propaga- 
tor. This property is to be contrasted with that of the standard or Symanzik 
improved actions where the spectrum is distorted by 0(a 2 ), or 0(a 4 ) cut-off 
effects, respectively. A distorted spectrum implies distorted thermodynamics. 
Large cut-off effects are seen in the free energy at high temperatures both in 
perturbation theory and in numerical simulations fl2] ] . 

The spectrum is determined by the singularities of the propagator D fll/ (q). 
Consider this propagator after n RG steps as given by Eq. (A7) in Appendix A. 
The ~ part of D^} (q) gives for large n the term 



2ttZ 



2" 



(24) 



The tensor QS n ^ is a smearing function which builds up the blocked link after 
n RG steps out of the original link variables. This is a local function in config- 
uration space whose Fourier transform is smooth, flffliq — 0) = <^„. Eq. ( p4| ) 
has a pole at fco = j fc] , where k = q + 2-kI. For any given q, qi £ (0,27r) 
the summation over I gives a tower of poles in D^^q) which reproduce the full 
relativistic spectrum of the continuum theory. 



The terms proportional to 1/k in Eq. (A7) are constructed out of Sl^y and 
do not change this conclusion^. Similarly, the part of D^ v proportional to the 
gauge fixing parameter a is removed at the end and does not influence the form 
of p FP (q) (cf. Appendix A). 

6 This part of the FP propagator is more complicated than the one in the non— linear a— 
model because with our block transformations a link variable contributes to several block 
averages. 
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In Fig. [j] we compare the spectrum of the quadratic parts of the Wilson, 
Symanzik improved and FP actions. Our parameterization of the Symanzik 
actions follows the notation of Weisz |?J . For the Wilson and Symanzik actions 
the momentum in Fig. [?] is restricted to lie along one of the lattice axes. The 
energy- momentum relation for the action labelled SI becomes complex at large 
k; the dotted line shows its real part. 




Figure 7: The spectrum of the quadratic part for the Wilson action (W), the 
Symanzik improved action for two different choices of the coefficients: SI: c\ = 
-1/12, c 2 +c 3 = 0, S2: ci = 0, c 2 + c 3 = 1/12, and the FP action. 
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6 The connecting tensor Z^ v {n) and the FP ac- 
tion at the cubic level 



In an SU(N) non Abelian theory gauge symmetry implies the presence of terms 
in the action which are cubic in the vector potentials. Gauge symmetry alone, 
however, does not fix the cubic part of the action. There are combinations of 
Wilson loops whose contribution to the quadratic part is zero while to the cubic 
part is not. The cubic part of the FP action is, however, uniquely given by the 
FP equation Eq. @. 

One expands Eq. ( |To| ) up to terms cubic in the coarse (B^) and fine (A^) 
vector potentials. The action is written in the form 

2N ■ S(U) =\Ys P^(r)A»(n + r)A» + (25) 

where in terms of the SU(N) generators, 

[T a ,T b ] = if abc T c , TrT a T b = ^6 ab . (26) 



The terms proportional to n in Eq. (|10|) are expanded to cubic order as well. 
The quadratic and cubic coupling constant tensors Pp, v (r) and c iivp {rx,r2) are 
determined then by the FP condition Eq. (|To|). 

After introducing a gauge fixing term in the action and solving the quadratic 
problem one obtains the minimizing fine vector potentials as a function of the 
coarse fields (Appendix A). In configuration space this relation can be written 
in the form 

A«(n) = Z ^ n - ^' B )B a p (n' B ). (27) 

Since the minimizing configuration is fixed only up to gauge transformations (cf. 
Section 2), the very definition of the relation in Eq. (|27]) requires gauge fixing. 
We used the type of gauge fixing given in Eq. ( |22| ) and defined the connecting 
tensor in the a — ► oo limit. Using eqs. (A5),(A6) one obtains 

1 1 1 

Z w {n - 2n' B ) = Jim ±- £ ^ £ e <(i*B-HrD(»-2»i) x (28) 

,0 ° B k B 1=0 



HP 



15 



n 


7 1 

z oo 


7 ii 

^00 


n 


7 1 


7 11 

^10 





0.4031 


0.4442 





0.0878 


0.0751 


10 


0.1756 


0.1413 


10 


-0.0408 


-0.0475 


11 


0.0610 


0.0623 


2 10 1 


0.0226 


0.0153 


-110 


0.0514 


0.0526 


3 10 


0.0217 


0.0148 


111 


0.0339 


0.0345 









Table 3: Some of the elements of the tensor Z pv which gives the minimizing 
configuration in terms of the coarse vector potential for the block transformation 
type I and type II. 



Note that sites on the fine lattice can be labeled as n = 2ns + m, m p = 0,1. In 
Fourier space one obtains 



The tensor Z pv is an important quantity which enters not only the cubic 
calculation, but also the construction of the fixed point operators, like the FP 
Polyakov loop (cf. section 8). The connecting function decays rapidly with 
distance; the minimizing configuration is determined locally by the coarse vector 
potential B p . In Table [|a few leading elements of Zoo and Z\q are collected for 
the two different block transformations. A permutation of the coordinates and 
the (somewhat unusual) reflection properties of Z pu (Appendix B) help one to 
find its value for further related fx, v and n values. 



The minimizing vector potential is given by Eq. (|27| ) only in the quadratic 
approximation. Nevertheless, this expression can also be used when we need the 
minimum in Eq. (^) up to cubic order. The error made is only quartic. This 
way all the cubic terms in Eq. (|ic|) will be expressed in terms of the coarse fields 
B® and one obtains the following recursion relation for the cubic coefficients 



(n' - n,n" - n)Z^^{n' - ri)Z v > v (n" - r 2 )Z plp (n) 



(29) 



n,n' ,n" 



where H pl/p comes from expanding Eq. (Jig) up to cubic order. The form of H t 



depends on the block transformation. We discuss here neither the explicit form 
of this tensor nor the details of solving Eq. ( |29| ) for the FP cubic couplings. 
In Table [| we give some of the largest elements of c pvp {r\,r2) for the block 
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n 






T2 




cooo(ri,r 2 ) 


-2 








-1 








-0.019 


-1 


-1 





-1 








-0.012 


-1 


-2 





-1 


-1 





-0.005 


-2 


-1 





-1 


-1 





0.002 




n 










cioo(ri,r 2 ) 





-1 








-1 





-0.598 





-1 -1 








-1 





0.027 





-1 





-1 


-1 





-0.027 





-2 








-2 





0.025 




n 






T2 




C2w{ri,r 2 ) 





-1 -1 








-1 





-0.062 





-1 -1 


-1 





-1 


-1 


-0.007 


-1 


-1 








-1 -1 





-0.004 





-1 -1 








-1 


-1 


-0.004 



Table 4: Some of the largest cubic couplings for the type I RG transformation. 

transformation of type I (for type II we have not yet studied the cubic problem). 
Other related elements can be obtained by using the symmetry properties of the 
cubic coupling tensor. 

We shall return to the cubic couplings in paper II when we discuss the general 
problem of the parametrization of the FP action in terms of loops. 

7 Cut— off independence of physical predictions 
of the FP action in 1— loop perturbation the- 
ory 

As discussed in Section 4, in leading order (tree level) perturbation theory the 
FP action gives cut-off independent physical predictions. Formal RG consid- 
erations suggest that the renormalized trajectory coincides with the FP action 
even at 1-loop level. That would imply that the predictions using the FP action 
are free of lattice artifacts even in 1-loop perturbation theory. Said differently: 
the predictions of the classical perfect action (FP action) contain no 0(a") nor 
O (g 2 a n In a, g 2 a n ) corrections for any n > 1. 

This is a strong statement which, unfortunately, we are not able to support 
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by a rigorous proof. We shall present the formal RG arguments for an infinite 
system which show that the only effect of the perturbative g 2 corrections is 
to make the marginal operator S FP (U) weakly relevant, i.e. the coupling /3 = 
2N/g 2 begins to move, but the form of the action on the renormalized trajectory 
remains unchanged in this order: 

f3S FP (U) — ► f3'S FP (V) , (30) 

under a RG transformation, to 1-loop level. Here /3' = /3 — A/3, where A/3 is 
fixed by the first universal coefficient of the /3-function. (For a scale 2 block 
transformation A/3 = 0.579 in SU(3).) An equation analogous to Eq. ( ^0|) 
occurred in ref. |jl3f where Wilson referred to formal RG arguments which were 
not, however, presented there. 



Stronger arguments come from an explicit perturbative calculation with the 
FP action of the non-linear c-model in d = 2 In ref. || the mass gap m(L) 
was calculated in a finite periodic box of size L. The prediction of the standard 
action is contaminated by power like cut-off effects (a 2 /L 2 ) n , n = 1, 2, .... These 
type of cut-off effects are eliminated when the FP action is used. More precisely, 
in the actual calculation they are reduced to the level of the numerical precision^. 
For example, in the range of L ~ 5 — 10 the cut-off effect is reduced by 5 orders 
of magnitude relative to that of the standard action. On the other hand, there 
exists a small non-power like cut-off effect, which is seen at L = 2, and goes 
to zero so rapidly that for L > 3 it disappears in the numerical errors. This 
behaviour can be associated with the small, but finite, range of the FP action. 

Let us present now the formal RG arguments. Consider a theory which has 
one marginal operator Ri(U) and many irrelevant operators R n (U), n = 2, . . .. 
The action is 

S = R 1 (U) + c 2 R 2 {U) + ... (31) 

and the Boltzmann factor is e~^ s ^ u \ (3 — 2N/g 2 for SU(N) gauge theories, 
while R\{U) is the FP action S FP (U) in our earlier notations. Under an RG 
transformation 

S — > S' = R!(V) + c' 2 R 2 (V) + . . . 

(3 — > 0'. (32) 

If the theory is asymptotically free the critical surface is at g = 0. A linearized 
scale s transformation at the fixed point is described by 

9^9 (33) 
c i Vj ( 34 ) 

7 The numerical errors are dominated by the errors in the couplings of the quartic part of 
the action. 
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where Xj — s kj < 1 since the Rj's are irrelevant. As the fixed point is at g = 0, 
the exponents kj are given by dimensional analysis. 

Now we consider general RG flows near the FP. For simplicity we introduce 
the notation c\ = g 2 . 

c[ = c^+a^c^+Oic 3 ) (35) 
4 = XkC k + 4fc iCj +0(c 3 ). (36) 

If we start with the action 

S FP = R U (37) 

Ci small, c 2 = c 3 = ... = 0, i. c. the classical fixed point action, then after one 
RG step 

c'j = a +a[ 1 1 ) c 2 1 + 0(c 3 ) (38) 

4 = a^cj + Oic 3 ), (39) 



or 

^ 2 S FP - - \ 1} ~ (S FP + a$g*R 2 + a$ g *R 3 + . . .) (40) 
9 9 + a n9 



Noticing that the prcfactor can be written 



> AiT7 = i-^ )+0(32) (41) 

9 2 + u\i9 A 9 

we see that to leading order in g 2 the effect of the RG transformation is only to 
shift the coupling 



PS FP -»(/3-aiV) [s FP + ^R 2 + ..^ , 



(42) 



and that the action itself is changed only in order g 4 . The shift in /?, is 
equal to A/3 introduced earlier. 



8 Fixed point operators 



The action on the renormalized trajectory is perfect — it has only physical 
excitations. In a Green's function only physical states come in as intermediate 
states independent of the operator used. In spite of that Green's functions (of 
fields, currents, etc.) show cut-off effects in general. Consider, for example a 
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free scalar field. As Eq. (21) of ref. 0] shows the perfect action on the lattice 
is described by a lattice field which is the integral of the continuum field over a 
hypercube. Therefore, the two-point function of the lattice field in the perfect 
action is equal to the two-point function of the continuum field averaged over 
the hypercubes around the the two lattice points in question. This averaging 
brings in a trivial rotation symmetry breaking. The two-point function of the 
field in the perfect action is not rotationally invariant, although only physical 
states propagate. 

A more interesting example is the static quark creation operator represented 
conventionally by a Polyakov loop (or by one of the lines of an elongated Wilson 
loop) The two-point function gives the quark-antiquark potential. The potential 
measured with the Wilson action at distances a , \[2a ,V3a, 2a, ... is non-smooth 
even at large beta values — it reflects the hypercubic lattice structure rather 
than physics. This is unfortunate since in numerical simulations these points lie 
in the physically interesting region. In addition, these are the points with the 
smallest statistical errors. One expects that simple Polyakov loops constructed 
in terms of the fields of the perfect action behave better, since the cut-off effects 
due to the distorted spectrum are eliminated. However, in order to get rid of 
the remaining cut-off effects we have to modify the operator itself. 

The first step in this direction is to construct FP operators along the lines 
by which the FP action was constructed. In Section 8.1 we use the example 
of the Polyakov loop to illustrate how to obtain the corresponding classical FP 
equation. In Section 8.2 we turn to the simple example of a scalar field theory 
in d = 2 and construct the FP field on the lattice. We show that the two-point 
function of the FP field is free of any cut-off effects of order (a/r) n . On the 
other hand, the FP field two-point function differs from the exact continuum 
propagator — (27r) lnr + const by terms which decrease exponentially with 
r. In Section 8.3 we solve the FP equation for the Polyakov loop in quadratic 
approximation in the vector potentials, calculate the corresponding qq potential 
(which is — (Airr)^ 1 in the continuum limit) and compare it with that of the 
standard and Symanzik improved actions. 



8.1 The FP Polyakov loop 

Consider the correlator (L(n)L^ (m))\ SFP , where L(n) is a Polyakov loop run- 
ning in the temporal direction at the spatial point n — (rii, ri2, n^) (with trace). 
Perform a RG transformation as in Eq. (Q) 

J DUe-P( sFP ^ +T ^)L(U; n)L\U; rh) . (43) 
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The configuration {V} on the coarse lattice is arbitrary; it can be smooth or 
rough. In the (3 — > oo limit (with {V} fixeoQ) Eq. ( p3| ) is a saddle point problem. 
One obtains 

e -0S^(v) Ll(y . HB)L n (v . Ab) 5 (44) 

where 

L'(V-n B ) = L{U{V)-n = 2n B ) . (45) 

Here U = U (V) is the minimizing configuration of Eq. (||) and we assume that 
all the coordinates of the two points n and m are even, n = 2ris, fn = 2ms. 
The FP Polyakov loop operator is defined as 

\L FP (v; n B = 1^ = L FP (U(V);n) , (46) 

for arbitrary configurations {V}. The eigenvalue A is expected to be 1 for the FP 
Polyakov loop. Indeed, the formal limit of the Polyakov loop is / d A x current x 
vector potential for a heavy particle, and this is a standard, dimensionlcss part 
of the action. This expectation is confirmed solving the FP equation Eq. ( (f6| ) 
at the quadratic level (Section 8.3). 



8.2 The fixed point field in a free scalar theory 

We shall construct the FP field operator in a d = 2 massless scalar field theory 
and investigate its properties. It is easy to handle this problem analytically 
and to demonstrate the essential features of the FP operator. Let us denote the 
scalar fields before and after the RG transformation by 7r„ and \n B > respectively. 
The RG transformation is defined as 

CXP ] 2 ^ P'^ B ^ n B^B+r B \ = (47) 
I n B ,r B ) 

C I Dvrexp I ~ ^ p(r)n n n n+r - 2k ( X n B - j ^ nn ) 

I n,r n B \ n£n B / 

where k is a free parameter. The corresponding FP action has been constructed 
by Wilson and Bell juj a long time ago. This free field problem is also relevant 
for the FP action of the non-linear cr-model [||. We shall follow the notations 
of ref. H . The FP solution of Eq. MTn in Fourier space has the form 



8 Observe that this limit does not correspond to the standard perturbative expansion where 
the relevant {V} configurations are smooth, 1-V ~ /J" 1 / 2 . 
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where the summation is over the integer vector I — (lo,li) and g, G (0, 2ir). A 
Gaussian integral is equivalent to a minimization problem. The FP equation 
for Eq. ( f47| ) can also be written as 



mm ■ 



^ Y P FP ( r B)Xn B Xn B +r B = (49) 
n,r n B \ n^riB / 



Let us denote the minimizing configuration by 7f„(x). One can repeat the min- 
imization procedure by going to even finer lattices, i.e. making a multigrid 
minimization, in principle to arbitrary depth. Then the minimizing field on the 
finest lattice at some point x B (measured in units of the coarse lattice) is given 
by a function &(\,x B ) of the coarse field. Clearly, it gives the solution to the 
fixed point problem: 

\<5> FP ( X ,x B ) = <$> FP {n( X )-2x B + 1 -), (50) 

where A = 1 as expected for a dimensionless scalar field in d = 2. The FP field 
operator Q FP (ir;x) is a linear combination of the fields 7T n /. The procedure 
of finding the FP field is simplified by the observation that k RG steps are 
equivalent to a single blocking step with a scale factor of L — 2 k with somewhat 
modified parameters p"4j : 

CXP ] ~\ Y P FP ( r B)Xn B Xn B +rB \ = (51) 

J L>7rcxp I - i^p FP (r)7r„7r„ +r [ Xn B ~ hk Y 71-71 ) f ' 

I n,r riB \ n£n B / \ 

where the summation n 6 n B goes over the L 2 = 4 fe fine points of the block 
labeled by n B and 

°* = tzi7Z5' hk = h L = 2k - (52) 

The minimizing configuration 7f n (x) of the quadratic form in the exponent of 
Eq. ( |5l|) gives then the FP operator in the k — > 00 limit 

® FP (x,n B ) = hm 7f n ( x ). (53) 

k—>oo 



c 
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Here we take0 m = L(n B i + |)- I 11 Fourier space we have (jl4| 

ff(g) = ±ut(q)^^ X (Lq) (54) 



1 - in, L 



where 

= II 7=3*- ( 55 ) 

Taking the fc — > oo limit one obtains 

$ FP (n B ) = ]T C(n B - n' B ) Xn > B , (56) 



or in Fourier space 

<t> FP (q B ) =C(q B )x(q B ), (57) 

where 

no 1 -r^r sin f^r + 7rZ) . 

1— — OQ 2 — V 2 



Let us demonstrate now that the propagator of the FP field & FP is free of 
those cut-off effects which go to zero as an inverse power of the distance. We 
show first that the propagator of the X ne ld itself has such cut-off effects in 
spite of the fact that the spectrum is exact. Using Eq. ( jig) and introducing 
k = q + 2irl, we get 



(XnXo) 



ikn 



d 2 k e 



-oo (27T) 2 k 2 



n 



• 2 kj 

sin -j- 

WW 



n ^ 0. 



(59) 



Consider \n\ 3> 1 (i.e. low momenta) and expand in k: 



(XnXo) 



d 2 k e lkn f 



-oo (27T) 2 k 2 



12 



— (h 2 ) 2 — 

360 v ' f440 



(k 4 



(60) 



The (k 2 ) s , s > 1 terms in the bracket do not give long range contributions. On 
the other hand, the rotationally non-invariant term kg + k\ (since it does not 
cancel the 1/fc 2 singularity) leads to cut-off dependent long range corrections 



1 UnU 



2^2 



0"1 



240tt (n 2 ) 2 ' 30tt (n 2 ) 4 



(61) 



9 For the given block transformation the sites in the block ub are given by n = Lub + m 
where m; = 0, 1, . . . , L — 1 for the components i = 0, 1 and the coarse field is coupled to the 



average over this region. Hence the choice nij 
of the averaging region. 



L/2 represents the site sitting in the middle 
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Consider the propagator of the FP field in Eq. ( p7| ) 

^ = Pfy- (62) 

For small q we have 



1 i JL S in 2 f 

fpT\ 2 1 1 / /o\2 + re S- terms ' 63 

p FP (q) 9 2 f = g (*/ 2 ) 2 

-i 2 



|C( g )| 2 ^(p^(g)) 2 
which shows that 



1 • Qi 

sm tt 



i=0 

1 



(ft/2) 



^ FP (<z)^^(l + 9 2 i?(9)) , (64) 

where is regular at q = 0. In D FP (q) the singular 1/g 2 • X)i(ft) fe i & > 1 
terms cancel, therefore no cut-off dependent corrections enter which vanish like 
an inverse power of n. 

The FP propagator has exponentially decaying cut-off dependent correc- 
tions, however. The following considerations not only illustrate this statement 
but give a general argument why the power corrections are canceled. 



Collect the terms quadratic in tt in Eq. (51) and introduce the notation 

Q(n, n')ir n TT n > = ^ p F p (r)TT n n n+r + a k b 2 k ^ 7r n J . (65) 

n,n' n,r n B \n£n B / 

The exponent in Eq. (^l|) then has the form 

^2Q(n,n')(-K n -Tt n (x))(lT n ' -7f„'(x))+ ^ P FP ( r B)Xn B Xn B +r 



(66) 

Let n — 2 k {iiB + \) and n' = 2 k (n' B + i). For given ns, n' B and k — > oo 



A cont O - n') = i / Dinr n Tr n , exp i -~ p FP (r)7r m 7r m+r 1 



(67) 



is the continuum propagator — (2-7r) 1 In \n\ + const. Insert 1 in the path integral 
in Eq. (|67]) using 

1 ~ I Dxcxp I -7; a k^2 ( Xrn B - b k ^ K rn) > (68) 
I m B \ mEm B / J 
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Using Eq. ( pq ) and introducing the new integration variables 

= 7T„ - Tt n {x) > (69) 

we get 

A^(n-n') = ($ FP (x,n B )0> FP (x,^)>| pFP + (VmMIq . (70) 

Note that A cont (n - n') = A cont (n B - n' B ) + const. The difference between 
the continuum and FP propagators is the inverse of Q defined in Eq. (fiBf). 
The tensor Q(n,n') connects fields over a distance ~ iafinc = &coaree so ^ is 
expected to produce correlations at least over such distances. On the other hand, 
comparing eqs. ( |65|) and ([5l|), {i> n i>n')\Q can be interpreted as propagation in 
the presence of the external fields Xm B which are put now to zero at every 
coarse lattice point tub- This forces the block averages J2 m em B n m to fluctuate 
around zero. Such an external field problem is expected to produce short ranged, 
exponentially decaying correlations^]. 



8.3 The FP Polyakov loop at the quadratic level and the 
qq potential 

We solve now the FP equation for the Polyakov loop Eq. ( ff6|) by iteration in 
quadratic approximation in the vector potentials. At the quadratic level we 
have 

L(U,n;k) = N-^ ^ W^{ft' - ft,n" -n,n' - n'^k)A^{n')A^(n") , (71) 

n' ,n" 

where k keeps track on the number of iterations performed, n — (no,n) , etc. 
and a summation over the Lorentz and colour indices is understood. In Eq. ( |7l|) 
time translation invariance is used. W is the unknown function to be deter- 
mined. Using the relation Eq. (|27]) between the coarse vector potential B and 
the minimizing A one obtains the following recursion relation for W: 

W pa (n' B ,n' B ,n' B0 -n' B0 ;k + l) = (72) 
W„ v (n', n", n' - <; k)Z w (ri - 2n' B )Z vrT (n" - 2n B ) . 

n' ,n" 

The starting value (corresponding to a Polyakov loop at ft = 0) is 

wv^g ) - { ; tt jzz "• s ' = s " = m 

10 Actually, this is the basic assumptions of the theory of RG transformations. 
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Due to the fact that W(., k = 0) is independent of n' — n' ', W(., k) is time 
independent, too. In addition, W(.,k = 0) is factorized in n' and ft", from 
which the same follows for arbitrary k. This is a great simplification. Using the 
reflection properties of Z (Appendix B) one obtains 

Z Ol (n) = 0, i= 1,2,3, (74) 

where 

Z^ v {n) = ^ z iiu{n, n ) . (75) 
no 



Eq. ( |7l| ) implies that only the fjb = v = component of W is non-zero for 
any k. One finds: 

L FP (U, 0) = N - I w(n) ^ Ag(n , n)] (76) 

plus terms higher order in A, where w(n) is the A = 1 eigenvector of the following 
eigenvalue equation: 

^2 Z o(n - 2n B )w(n) = Xw(n B ) ■ (77) 

n 



One can solve the eigenvalue equation by iteration starting with w(n) = 5^ - 
The iteration will automatically project out the eigenvector with the largest 
eigenvalue A = 1. Given Z, this is a trivial problem. The error in w(n) is 
dominated by the error in the connecting function Z. For the case of type I RG 
transformation the result is shown in Table 5. The other elements are smaller 
than 10~ 5 . The eigenvector w(n) in case of type II RG transformation looks 
rather similar |jl^|. It follows from Eq. j77| ) and Eq. (B2) in Appendix B that 
w(n) is reflection symmetric under rij — > — n,, i = 1,2,3. The eigenvalue A is 
unity in the quadratic approximation, as expected. 



Let us consider now the correlation function of two FP Polyakov loops 

(L FP (U;n)L FP (U;0)) (78) 

expanding L FP up to quadratic order in the vector potentials. The connected 
part is proportional to the square of the lowest order quark-antiquark potential 
[pi. Using Eq. @ we get 



V(r) = 



2 71 



d 3 q 
(2nf 



"w(q)w(-q)D QQ (qo = 0,q), 



(79) 
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n 


w(n) 




71 


w(n) 







0.8951 




12 


0.1312 x 10" 


-3 


1 1 


0.5626 x 10" 


-2 


2 2 


0.3763 x 10 


-4 


111 


0.1729 x 10" 


-2 


1 3 


0.2709 x 10 


-4 


1 


0.1552 x 10" 


-2 


1 2 2 


0.1806 x 10 


-4 


2 


0.9364 x 10" 


-3 


1 1 3 


0.1124 x 10 


-4 


1 1 2 


0.1425 x 10" 


-3 


2 3 


0.3954 x 10" 


-5 



Table 5: The weights w(n) for the FP Polyakov loop. 

where D is the propagator of the FP quadratic action (Section 4 and Appendix 
A) and 

w (q) = J? e -^ a w(n) (80) 

n 

is the Fourier transform of the FP Polyakov loop. In the continuum limit we 
should get V{r) = — (47ir) — 1 . 



We demonstrate now that V(n) in Eq. ([79j) has no cut-off effects which 
for large \n\ decrease as an inverse power of \ft\. Said differently, the potential 
derived from the FP Polyakov loop is a solution of the (on-mass-shell) Symanzik 
program through all orders in a 2 . The argument follows the same lines as in a 
free scalar theory, eqs. ( |6^ - |64[ ). Separate first the singular part of D 00 (q =0,q) 
at small q. As Eq. (A7) in Appendix A shows, the singular contribution to 
after n iterations comes from the I = term 



D$(q o = 0,q)= n$(2- n q) 



>(«) 



1 



regular terms, 



(81) 



where we used E q. (|22|) and the fact that only the 00 component of fi survives 
at qo = as eqs. (17, 1J), and (A8) show. The presence of the non-rotationally 
invariant residue in Eq. ( pi] ) is responsible for the remaining cut-off effects 
decaying as a power of \n\ if we use simple Polyakov loops (w(q) = 1). This 
dangerous residue is cancelled by w(q)w(—q) in the FP loop correlators, however. 



Using Eq. (|77|) we get the recursion relation for the eigenvector with A = 1 
(the argument go = is not indicated explicitly in the following equations): 



h=0 



2irl 



,(n) 



2ttZ 



(82) 
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where Zoo is obtained from Eq. 



,(n) ( q + 2nl \ _ (n) ( q + 2irl 



J Q0 



q + 2nl 



-on 



Di n o +1 \q) 



(83) 



Iterating Eq. ( p2\ ) leads to 



- (n+1) ® = iK +1) ®l _i E fi 



2"-l 



)(«)* 
'00 



g + 27r/ \ ^(o) ( q + 2irl 



-'no 



z,=o 



For large n eqs. (gj),(|8lj) and (gg) give 



w {n+1 \q)=f 



1 



2"-l 



(n)* 
00 



2" 



(84) 



<f + 2-nl \ 1 



;=o 



(<?+ 2tt0 2 



(85) 

Here and below we denote by R(q) a generic term regular at q = 0. Only the 
? = term gives relevant contribution 



W {n+1) {q) = 



1 



(2"»g) 



a +g«i2(g)| (fi^dj+M?)] . (86) 



Combining eqs. (|79|), ( pl| ) and ( |86| ) we get 

f2T A .,.,v • 



V(f) 



which we wanted to show. 



o 



(2tt) 3 <f 



i + ^fo)) , 



(87) 



The considerations above say nothing about the exponentially decaying cut- 
off effects. For a graphical illustration in Figs. [j], || and || we plot the perturbative 
potential rV{r) versus r for the FP, Symanzik and Wilson actions. 



9 Outlook 



The most relevant task for future progress is to construct the FP action for full 
QCD including fermions. The free fermion case and the anomaly problem in 
the Schwinger model have been considered in this context Jl7[ . The QCD FP 
action is quadratic in the fermion fields, which is an important simplification 
and should make this problem feasible. 
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A Appendix 



In this Appendix we discuss the steps of solving Eq. (|T^) for the FP p^f . 
In the first term on the r.h.s. of Eq. ( |T9| ) we write 



V B 16 • 

k B I 



where = 0,1. Taking the derivative of the r.h.s. with respect to A v (—q) for 
q = %y + ttI one obtains 



a v + in) + K [d + ^ w t + ^] ^ (r y to ) - bm) = o. 

^ (A2) 

It is assumed that some gauge fixing term is introduced and so p has an inverse 
which is denoted here by D. First one solves Eq. (A2) for T^. After multiplying 
with w lv + 7rZ) and summing over v and / we get 



M{q B ) + ~ 

K 



Mp, p {q B )B p (q B ) 



7i» 



where 



(A3) 



(A4) 



After substituting Eq. (A3) into Eq. (A2), the minimizing is obtained as a 
function of the coarse field B v 



l 



pp \K 



+ M{q B )\ B u {q B ). 



pv 



{AS) 
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We can write now Eq. (A5) into Eq. (|T^) and find after some algebra 

d '^Mb) = M^{q B ) + -<W 



(A6) 



which is the result quoted in Eq. (E 



In iterating Eq. ( pp| ) let us start with some D^J. We might take the form 
given in Eq. d23). After n RG steps we get 



where 



and 



q + 2ixl 



D (0) 



q + 2nl \ t ( q + 2irl 



fi$ (2" J 'fc) = (w(2- 1 fe)a;(2- 2 fc) . . . «(2"^)) 



//// 



^ n— 1 



(A7) 
(A8) 



2 J -1 



Q$(9) = iJ E [fi W (2^(g + 27rO)n Wt (2- J ' (g + 2ttZ)) 



(A9) 



One can prove eqs. (A7-A9) by induction. We introduced the following 
shorthand notations above 



2"-l 2-1 2-1 

E = E-E 

1=0 l =0 l 3 =0 



(AW) 



The tensor 0$, (k) in Eq. (A8) can be written in a form which is well suited 
for a numerical procedure. For both types of RG transformations Wpi,(k) in 
eqs. ( |l7|Jl^ ) can be written as 



^(fc) = ^a( 1 )(fe)V + ^6«(fc). 



One obtains then 



(All) 



n$(2-j*) = 2-'". 



, '_ = #_ (J')(2-J'jfe) + 2^'^^6W(2- J 'fc) . 
2-J'feu 2-Jfc„ 



where a^-* and satisfy the simple recursive relation 



0412) 
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b^(2- j k) = b^(2- 1 k)a iJ - 1) (2- j k) + (2~ j k) . (AL3) 

The equations above allow a sufficiently precise numerical determination of the 
fixed point propagator D^(q). 

Let us now discuss briefly the steps of removing the gauge fixing. The brute 
force method is to take the a — > oo limit of p^f = (D fp J ^ numerically. 
More elegant is to observe that the longitudinal part of the propagator remains 
longitudinal after RG transformations: 

D$ = £> (0) (<?)<W + «?^/ (0) (9) — D<$ = G$(g) + aq&frHq) ■ (AU) 

where G0 \q) is the part evolved from D^ Q \q)8 llv . Eq. (A14) can be shown 
easily by using the gauge symmetry relation satisfied by the tensor 

Using Eq. (A14) one obtains for large a 

/ [iv * a \ a / 

where p^J (q) is the gauge invariant part of the quadratic blocked action we are 
looking for. Introducing the projector 

V(<?) = <W-||| {An) 

it is easy to show that 

p fiV (q) = (T( q )G(q;r l )- 1 T(q))^ , (A18) 

where ^ 

G(q; v)„ = (T(q)G(q)T(q))^ + „Mi . (A19) 

Here 77 is an arbitrary non-zero parameter and p in Eq. (A18) is independent 

of 7]. 



B Appendix 

Using the definitions Eq. ( |l3| ) and Eq. ( p7| ) for p^ v {r) and Z M „(r), respectively, 
it is easy to derive the following transformation properties under reflection of 
the coordinates. 



31 



Consider the index combinations [iv = 00 and 10. The other cases can be 
obtained by permuting the coordinates. We get (r = (ro, ri, r 2 , r 3 )) 

Poo(ro,ri,r 2 ,r 3 ) = poo(±r , ±n, ±r 2 , ±r 3 ) , 

Pio(t"o, ri, r 2 , r 3 ) = -pio(-r-o + 1, n, ±r 2 , ±r 3 ) , (Bl) 
pio^o,^^^) = -pio(r-o,-r-i - l,±r 2 ,±r 3 ), 

Z 00 (r ,ri, r 2 ,r 3 ) = Z 00 (r , ±n, ±r 2 , ±r 3 ) , 
^oo(r"o,ri, r 2 ,r 3 ) = Z 00 (-r + 1, ±n, ±r 2 , ±r 3 ) , 
^io(r-o,n,r- 2 ,r 3 ) = -Z 10 (-r + 2,n,±r 2 ,±r 3 ) , (52) 
^io(r-o,n,r- 2 ,r 3 ) = -^i (r ,-ri - 1, ±r 2 , ±r 3 ) . 
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